Magnetic confinement of massless Dirac fermions in graphene 



o 
o 

c 



A. Dc Martino, L. Dcll'Anna, and R. Egger 
Institut fur Theoretische Physik, Heinrich-Heine-Universitat, D-40225 Diisseldorf, Germany 

(Dated: February 4, 2008) 

Due to Klein tunneling, electrostatic potentials are unable to confine Dirac electrons. We show 
that it is possible to confine massless Dirac fermions in a monolayer graphene sheet by inhomoge- 
neous magnetic fields. This allows one to design mesoscopic structures in graphene by magnetic 
barriers, e.g. quantum dots or quantum point contacts. 

PACS numbers: 73.21.-b, 73.63.-b, 75.70.Ak 



O 



-a 
o 

(N 
> 
O 
On 
(N 
O 

o 

-I— > 
c3 



i 

C 

O 

o 



X 



The successful preparation of monolayer graphene 
films 0, 0, [|| has recently generated a lot of excite- 
ment and allows one to directly probe the physics of 
two-dimensional (2D) Dirac- Weyl fermions. The mass- 
less Dirac spectrum at low energy scales is caused by 
the sublattice structure (the basis of graphene's honey- 
comb lattice contains two carbon atoms, giving rise to an 
isospin degree of freedom) combined with a special band 
structure, and has been verified experimentally P, [H, El- 
Besides the fundamental interest, graphene has also been 
suggested as a building block for future nanoelectronic 
devices Q. However, there is an interesting twist at that 
point, since Dirac fermions cannot immediately be con- 
fined by electrostatic potentials. In marked contrast to 
the Schrddinger case, Dirac fermions can penetrate high 
and wide electrostatic barriers with high transmission 
probability, in particular for normal incidence. This is 
often referred to as Klein tunneling [f| and can be under- 
stood by noting that under the barrier, the whole spec- 
trum is shifted upwards. Incoming electron-like quasi- 
particles can then efficiently tunnel through the barrier 
via empty states in the hole band, which are always avail- 
able since the Dirac spectrum is unbounded. In the con- 
text of carbon nanotubes and graphene, this effect was 
theoretically studied in Rcfs. [g, S Hi- The creation of 
useful mesoscopic structures, e.g. quantum dots or quan- 
tum point contacts, thus seems to encounter a major and 
fundamental obstacle, seriously limiting graphene's po- 
tential for applications. 

An obvious but rather crude way out of this dilemma is 
to mechanically cut samples into the desired shape. Al- 
ternatively, one could attempt to exploit the fact that 
suitable transverse states in a graphene strip may al- 
low one to circumvent Klein tunneling Q. Here we de- 
scribe a completely different and hitherto unnoticed way 
of confining Dirac- Weyl quasiparticles in graphene by 
magnetic barriers. Employing existing technology, the 
required inhomogencous static magnetic field configura- 
tions can be created using ferromagnetic layers located 
beneath the substrate on which the graphene layer is de- 
posited; for other possibilities, see Ref. [io| . Mesoscopic 
transport with magnetic barriers has been experimentally 
studied for the Schrddinger fermions realized in conven- 
tional semiconductor heterostructures, e.g. transport in 
the presence of magnetic barriers [ill and superlattices 
fl2| . magnetic edge states close to a magnetic step [Hj], 



and magnetically confined quantum dots or antidots [14| . 
Correspondingly, apart from one study of magnetic edge 
states in narrow-gap semiconductors [15| , model calcula- 
tions have only been carried out for Schrddinger fermions 
[ToL [l|| [13, EH; Here we formulate the theory of mag- 
netic barriers and magnetic quantum dots for the mass- 
less Dirac- Weyl fermions in graphene. With minor modi- 
fications, the theory also covers narrow-gap semiconduc- 
tors. 

We focus on a static orbital magnetic field ori- 
ented perpendicular to the graphene (x — y) plane, B = 
B(x,y)e z , and work on the simplest possible theory 
level (no disorder, no interactions), where the electronic 
spin degree of freedom can be disregarded. Moreover, 
we consider the physically relevant case of slow B(x,y) 
variations on the scale of the graphene lattice spacing 
(a = 0.246 nm). On low energy scales, a weakly doped 
(or undoped) graphene layer is described by two identical 
copies of the Dirac Hamiltonian, which remain decoupled 
in the presence of smoothly varying magnetic fields [2fJ. 
These two copies describe low-energy envelope states in 
the k • p approach [2l[ close to the two relevant K points 
in the hexagonal first Brillouin zone of graphene. For 
slowly varying B = rot A, we therefore need to study 
just one K point. The time-independent Dirac equation 
for the spinor ip(x,y) = (ip+,ip-) T at energy E = Vpe 
then reads (we put H = 1) 

(p+ ^ A (^y)) i>(x,y) = eip(x,y), (1) 

where the Fermi velocity is vp « 8 x 10 5 m/sec, the 
momentum operator is p = —i(d x ,d y ) T , and the 2x2 
Pauli matrices in a = (a x ,a y ) act in isospin space. The 
velocity operator follows from the Heiscnberg equation 
as v = vplj. In this paper, we discuss two prime ex- 
amples of interest based on the Dirac- Weyl Hamiltonian 
in Eq. ([1]), namely, (i) the magnetic barrier and (ii) a 
circularly symmetric magnetic quantum dot. 

For a magnetic barrier, the relevant physics is de- 
scribed by a magnetic field translationally invariant along 
the (say) y-direction, B(x,y) = B{x). Choosing the 
vector potential in the gauge A(x, y) = A(x)e y with 
d x A(x) — B(x), transverse momentum p y is conserved, 
and for given p y , Eq. (JTJ) leads to the coupled equations 



[d x ±p y ± (e/c)A(x)]ip±(x) = ieip T (x). 



(2) 
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These equations imply the decoupled ID 'Schrodingcr 
equations 



(3) 



with the Py-dependent effective potentials 

V±(x) = ±(e/c)d x A(x) + [ Py + (e/c)A(x)] 2 . (4) 

Let us then describe the solution of Eq. ([3|) for a square- 
well magnetic barrier, where B = B§e z (with constant 
Bo) within the strip —d < x < d but B = otherwise, 



B(x,y)=B 9{d 2 -x 2 ), 



(5) 



with the Heaviside step function 8. The sharp-edge form 
([5]) is appropriate when the Fermi wavelength Xp is para- 
metrically larger than the edge smearing length A s , while 
A s > a to ensure smoothness of B; otherwise scattering 
between the two K points takes place (23|. Note that 
here Xf ~ l/N is determined by the (inverse) Fermi mo- 
mentum of the Dirac quasiparticles, which is measured 
relative to the relevant K point, and the large momentum 
scale associated with the K point itself drops out com- 
pletely. With the magnetic length l B = ^c/eB , the 
vector potential is written as 



A(x) = -jtX 
ei B 




(0) 



Consider now an electron-like scattering state (e > 0) 
entering from the left side, with incoming momentum 
P = {Px,Py)- The incoming wave function is, up to an 
overall normalization, 



i>in(x) = [ Px +i( Pv -d/i%) J e ipxX , 
V IpI / 

where the shift in p y is due to our gauge choice for the 
vector potential. It is then convenient to parametrize the 
momenta as 



p x — e cos (/>, p y = e sin <ft + d/l\ 



(7) 



The gauge- invariant velocity is v = up (cos (f>, sin <fi) T , and 
therefore <j) is the kinematic incidence angle. The emer- 
gence angle <j>' at the right barrier, p' x = ecos</>', is ob- 
tained by exploiting conservation of p y , 



2d 



sin i 



sin < 



(8) 



Up to an overall normalization factor, the scattering state 
in the three regions is as follows. For x < — d, 



ip T (x) 



1 

oi4> 



1 



(9) 



with ^-dependent reflection amplitude r. In the bar- 
rier region \x\ < d, the solution is expressed in terms 
of parabolic cylinder functions D v [24| , 



D (el B )y2-l {±V2(x/l B +p y l B )) 



^D,.,_„ ,„ (±V2(x/l B + Py l B )) 



(10) 

with complex coefficients c±. Finally, for x > d, the 
transmitted wave is 

i>m(x) = t^pM ( e i' ) e ip '°> X (11) 

with transmission amplitude t. The transmission proba- 
bility T = \t\ 2 is then related to the reflection probability 
R = \r\ 2 by T + R = 1. Note that Eq. © implies that for 
certain incidence angles <f>, no transmission is possible. In 
fact, under the condition 



el B <d/l B , 



(12) 



every incoming state is reflected, regardless of the inci- 
dence angle <j>. In essence, all states with cyclotron radius 
(defined under the magnetic barrier) less than d will bend 
and exit backwards again. This illustrates our main find- 
ing: in contrast to electrostatic barriers, magnetic bar- 
riers are able to confine Dirac-Weyl quasiparticles. For 
sufficiently large barrier width 2d and/or field Bq, all rel- 
evant states will be reflected. Our analysis also shows 
that this conclusion is generic and does not depend on 
the particular choice (O for the barrier. 

If the condition (fT2"|) is not obeyed, the transmission 
probability T does not vanish in general. Its value follows 
by enforcing continuity of the wavefunction at x = ±<f. 
The solution of the resulting linear algebra problem yields 
the transmission amplitude in closed form, 



t = 



gi(px+pi)d J) 

V = (d B ) 2 e < ^'-^(uf«2 -<4%) 
— 2(vj t ~t> 2 ~ — ) + iV2el B 

where we use the shorthand notation 



uf = D 



(^b) 2 /2- 



_j [±V2(-d/l B +Pyl B )) , 



u± = D 



(eifl) 2 /2 



±V2(-d/l B +p y l B ) 



The related symbols uf,vf follow by letting — d — > d. 
The resulting transmission probability T{4>) = \t\ 2 is 
shown in Figurc[T]for several parameter values (e, d) out- 
side the perfectly reflecting regime specified in Eq. (Tj"2"j) . 
For a typical value of Bq = 4 T, the magnetic length is 
l B = 13 nm, and el B = 1 corresponds to E — 44 meV. 
Standard fabrication and doping techniques should thus 
be sufficient to enter the perfect reflection regime. On the 
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FIG. 1: (Color online) Polar graphs depicting the transmis- 
sion probability T(</>) for a magnetic barrier of width 2d at 
various energies e. The outermost semicircle corresponds to 
T — 1, the center to T = 0, with grid spacing 0.2. Angles 
between —tt/2 and +n/2 are shown, the angular grid spacing 
is 7r/6. (a) T as function of barrier width for fixed energy, 
eZs = 3.7. For cI/Ib > 3.7, the transmission is zero for all (f>. 
(b) Same as function of e for d/ls — 1-5. The transmission 
vanishes for ds < 1.5. 



other hand, for sufficiently high energy and/or narrow 
barriers, the ^-dependent transmission profile in Fig. [1] 
should be observable. 

We now turn to a discussion of a circularly symmet- 
ric magnetic quantum dot, defined by a radially inho- 
mogencous field B = B(r)e z . It is convenient to use 
complex variables z = x + iy (not to be confused with 
the z-direction) and z = x — iy, and the corresponding 
derivatives d = \{d x — id y ) and d = \{d x + id y ), where 
r = \fz~z. Writing in a similar manner A ~ A x 



-iA y and 



A = A x — iA y , the magnetic field is B(r) = —i(dA-dA), 
where (e/c)A = iip(r)/z and (e/c)A = —iip{r)/z. This 
gauge expresses the vector potential in terms of the mag- 
netic flux if(r) through a disc of radius r in units of the 
flux quantum hc/e, 



dr'r'B(r'). 



(14) 



Next, we recall that in 2D the group of rotations is 
50(2) ~ U(l), whose generator is the orbital angular 
momentum operator L = zd — zd, with eigenfunctions 
~ z m (integer m). For an isotropic field B(r), the opera- 
tor J = L + cr z /2 is conserved, i.e. eigenstates of Eq. (JXJ) 
are classified by the half- integer eigenvalue j = m ± 1/2 
of J, 



'</•+ 



1 



4> m {r)(z/r) m 
X m (r)(z/r) m+1 



(15) 



The Dirac equation ((T|) then reduces to a pair of radial 
ID equations for <p m (r) and Xm(f) (where r > and 
f' = df/dr), 



ip(r 



+ 1 + <p{r) 



(16) 



implying a second-order equation for the upper compo- 
nent of Eq. (TT5)) . 



-B(r) 



An = 0, (17) 



plus a similar equation for Xm(r). For e ^ 0, Xm directly 
follows from (f> m via Eq. (fTTJ)) . 

We now analyze a simple model for a magnetic quan- 
tum dot, where B = Bq outside a disk of radius R 
and zero inside, as previously considered for Schrodingcr 
fermions in Ref. 0. The flux dHJ) is with l B = y/c/eB 
given by 



R 2 



2l 2 



-R). 



(18) 



With normalization constant J\f„ 

6(r~R)R 2 /2l 2 B 



and m < 0, the states 







-V>(r)/2 





{z/r 2 ) 7 



represent zero-energy solutions of Eq. ([T]). The remain- 
ing eigenspectrum comes in pairs ±e, and we focus on 
the e > sector. Up to an overall normalization factor, 
Eq. (fl~7|) implies Bcsscl function solutions inside the dot, 
= J m for r < R. The general solution <^> > outside 
the dot (r > R) involves the degenerate hypergeometric 



functions $ and * [2J. With £ 



-s, 



where S = R 2 /2l 2 B is the missing flux through the dot, 
we obtain 



£|m|/2 e -£/2 ( ai $( Q!) i + |^| ; Q + a 2 *(a, 1 



N;£)) 

(19) 
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FIG. 2: (Color online) Low-energy eigenenergies (labeled by 
m) for a disk-like magnetic quantum dot in graphene versus 
missing flux S = R 2 /2l%. 



Here a\ t 2 are arbitrary complex coefficients, and energy 
is parameterized by 



l + m0(m) - (el B ) 2 /2. 



(20) 



Continuity of ip(r) at r = R now implies continuity of 
both 4> m (r) and 4>' m (r), see Eq. (fTtJ|) . The resulting two 
matching conditions then determine the possible eigen- 
states. 

Note that the well-known relativistic Landau levels (for 
R = 0) correspond to a = —n (with n = 0, 1,2, . . .) [2(j |. 
where $ and \I/ reduce to Laguerre polynomials. For 



finite R, the matching problem does not admit solutions 
with a = —n, and we thus consider a =/= —n. However, 
for a —n, $ has the asymptotic behavior $ ~ at 
£ — > oo, i.e. normalizability of ip necessarily requires a\ = 
in Eq. (|T9|) . One of the two conditions then fixes 02, 
and the other determines the quantization condition on 
the energy, 



1 - \ fh\9(-m)/S 
d 



d B J m+ i{el B V25) 



25 J m {d B V2S) 
ln*(a, 1 + |m|;£ = 5). 



(21) 



The numerical solution of Eq. (|2Tj) is possible using stan- 
dard root finding methods (bracketing and bisection) . In 
Figure O we show the solutions to Eq. (|21| with e > 
but below the lowest positive-energy bulk Landau level 
located at ds = \/2. Within the shown 5 range, for 
m 7^ 0, there is at most one solution with < el B < s/2, 
while for m = 0, we obtain two such solutions for 5 > 4. 
Depending on the missing flux S ~ R Bo, the energy 
levels of this 'Dirac dot' can be tuned almost at will. 

To conclude, we have described a new way of confin- 
ing Dirac- Weyl quasiparticles in graphene. We hope that 
our work will guide experimental efforts to the develop- 
ment of mesoscopic structures based on this novel mate- 
rial, and stimulate more theoretical work on the effects 
of magnetic barriers on Dirac fcrmions. 
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